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Sensitivity  Tests  of  a  Surface-Layer  Windflow 
Model  to  Effects  of  Stability  and  Vegetation 


I.  INTRODUCTION 

In  early  1984,  we  undertook  an  effort  at  AKGI.  to  acquire  an  appropriate  two- 

dimensional  windflow  model  for  use  in  conjunction  with  a  Gaussian  puff  diffusion 

1  2  .1 

model  developed  by  Fleischer  and  evaluated  and  modified  by  Kunkel.  ’ '  It  is  our 

intent  that  this  combined  prediction  system  for  toxic  chemical  transport  and  diffu¬ 
sion  be  adapted  for  use  by  the  Air  Weather  Service  (AWS)  over  regions  of  complex 
terrain. 

We  conducted  the  search  for  an  appropriate  windflow  model  with  the  following 
constraints  in  mind: 

(1)  The  model  must  tie  applied  over  complex  terrain. 


(Received  for  publication  16  October  1985) 

1.  Fleischer,  M.T.  (1980)  Si’ll. I.S  -  An  Kvaporation/Air  Dispersion  Mode)  for 

Chemical  Spills  on  l.and.  Shell  Development  Company,  W'esthollow  Research 
Center,  Houston,  Tex. 

2.  Kunkel,  H.A.  (1983)  Comparison  of  Evaporative  Source  Strength  Models  for 

Toxic  Chemical  Spills,'  AFGI.-fR-83-0.307,  AD  A1394.il. 

3.  Kunkel,  H.A.  (1984)  An  Evaluation  of  the  Ocean  Breezr/Drv  Gulch  Dispersion 

Model.  AFG I.-TR  -84-0313.  ADA  157165. 


2.  The  model  must  be  adaptable  to  either  a  small  minic  omputer  or  a  micro¬ 
computer. 

3.  The  model  must  be  adaptable  to  run  as  a  combined  system  with  the  AFCil. 
dispersion  model. 

4.  The  model  must  lie  adaptable  to  small  horizontal-grid  s pacings  (on  the  or¬ 
der  of  100-400  m)  on  a  horizontal  domain  approximate! v  10  by  ]0km  or  smaller. 

Given  these  constraints,  we  reduced  the  number  of  models  available  for  con¬ 
sideration  to  a  list  of  six.  A  summary  of  the  major  features  * 4 * 6  7 8 * 10 11 12’  *“ 

8  9 

and  the  disadvantages  *  of  each  model  considered  is  given  in  Table  1.  From  the 
information  contained  in  Table  1  and  the  prescribed  requirements  listed  above,  we 
selected  the  1  ,S.  Army  Atmospheric  Sciences  Laboratory  ( A S I . )  model  for  devel¬ 
opment  and  incorporation  into  the  combined  transport  and  diffusion  prediction  sys- 
t  e  m . 

This  report  summarizes  the  results  of  the  first  phase  of  the  windflow  model 
evaluation.  In  this  report,  we  will  discuss  the  modifications  made  to  the  ASI. 
model  to  produce  more  realistic  simulations  of  surface-layer  windflow.  We  will 
then  describe  the  sensitivity  tests  performed  with  the  model  using  simulated  mete¬ 
orological  data  and  detailed  vegetation  and  terrain  data  supplied  to  us  by  ASI.  for 
Fort  Folk,  La.  Model  sensitivity  was  tested  for  different  stability  conditions  (for 


4.  Hall,  J.  A.,  and  Johnson,  S,  A.  (1978)  Physically  Based  High  Resolution  Sur¬ 
face  Wind  and  Temperature  Analysi s  for  F  PAMS ,  \S I . - (  R  -TT- 004.3- 17 
F.S.  Army  Atmospheric  Sciences  Laboratory,  White  Sands  Missile  Range, 
N.  Mex. 

j.  Holdeman,  J.T.  (1982)  Modeling  the  tropospheric  wind  field  over  complex 
terrain.  Symposium  on  the  Composition  of  the  Non  Lrban  Troposphere, 
Williamsburg,  Va. 

6.  Dickerson,  M.  H.  (1978)  M  A  SCON  -  A  mass -consistent  atmospheric  flux 

model  for  regions  with  complex  terrain,  J.  Appl.  Mcteorol.  1 7:241  -253. 

7.  Dempsey,  I).  P. ,  and  Mass,  C.  F.  (11*83)  A  one-level  mesoscale  model  for 

complex  terrain,  Proc.  Sixth  Conf.  on  Numerical  Weather  Prediction,  A m ■ 
Meteorol.  Soc.  .  Boston,  Mass.  343-347. 

8.  Ilaltiner,  G.J.  (1971)  Numerical  Weather  Prediction.  Wiley,  New  York. 

0.  Walmsley,  J.  L.  (1983)  Modeling  of  Boundary  - 1  .aver  Wind  Flow  Above  Com- 
plex  Terrain:  A  Review,  Report  ARQH-35-flOa-T,  ATmospTierTc  ITriviron- 
rnent  Service,  Downsview,  Ontario. 

10.  Lavoie,  R.L.  (1972)  A  mesoscale  model  of  lake-effect  storms,  J.  \tmos. 

■Sri,  29:1025-1040. 

11.  Lavoie,  R.L.  (1974)  A  numerical  model  of  trade  wind  weather  on  Oahu,  Mon. 

Wea,  Rev.  1012:630-637. 

12.  Fox,  I).  G.,  Fosberg,  M.A.,  Marlatt,  W.K.,  and  Reeser,  VY.  (1976)  Analy¬ 

sis  of  mountain  air  quality,  Proc.  Third  Symposium  on  Atmospheric 
Turbulence,  Diffusion,  and  Air  Quality.  Am.  Mcteorol.  Soc.,  Boston, 
Mass.  470-473. 


example,  stable,  unstable,  neutral),  effects  of  horizontally  varying  roughness, 
and  reaction  to  the  presence  of  nonhomogeneous  vegetation  over  the  domain. 


2.  THE  WINDFLOW  MODEL 

In  this  section,  we  give  a  brief  description  of  the  AS1.  windflow  model  and  de¬ 
scribe  the  modifications  made  to  it  at  AFGL.  Details  of  the  model  are  given  by 

Hall  and  Johnson, A  Ohmstede,  ^  Ohmstede  and  Stenmark,  ^  Cionco,  ^  and 
1  6 

Amlicke  and  C  oleman. 

2. 1  General  Characteristics  of  the  Model 

Tlie  model  produces  a  two-dimensional  (x-y  plane)  surface-layer  windflow 

analysis  by  employing  Gauss's  Principle  of  Least  Constraints  to  adjust  an  initially 

uniform  windfield  to  conform  with  topography,  mass  conservation,  and  buoyancy 

forces.  Since  an  excellent  discussion  of  Gauss's  Principle  of  Least  C’onstraints 
1 7 

appears  in  I.anczos,  we  just  present  the  basic  model  constraint  equation: 

^volume  t  A  -  H1 J2  dV  *  aiinimum'  (1> 

where  A  includes  only  the  advection  terms  of  (lie  acceleration  and  B^  denotes  the 
surface  parallel  component  of  the  buoyancy  [(see  Kq.  (7)].  This  integral  is  taken 
over  the  volume  of  the  model  surface  layer,  which  is  divided  into  "flux  boxes”  of 
surface  normal  thickness  ,  conforming  to  the  warped  terrain  surfaces.  The 
flux  form  of  the  advection  in  Kq.  (1)  is  used,  and  a  terrain-following  coordinate 
system  is  used  to  include  effects  of  terrain  slopes  in  the  calculations.  The  total 
constraint  integral  is  expressed  as  a  sum  over  all  the  boxes  in  the  model  layer: 

1  1.  Ohmstede,  W.l).  ( 1 D8 1 )  The  parameterization  of  battlefield  dispersion-new 
frontiers.  Workshop  on  the  Parameterization  of  Mixed-Layer  Diffusion, 

Las  Cruces,  Mex.  279-287. 

14.  Ohmstede,  W  .  B. ,  and  Stenmark,  K.  B,  (19111)  A  model  for  characterizing 
transport  and  diffusion  of  air  pollution  in  the  battlefield  environment. 
Workshop  on  the  Parameterization  of  Mixed- 1  .aver  Diffusion,  Las  Cruces, 
N.  Mex.  416-422. 

la.  Cionco,  R.  M.  (19111)  A  meteorological  approach  to  chemical  defense  over 
complex  terrain  with  vegetation.  Workshop  on  the  Parameterization  of 
Mixed-Layer  Diffusion,  Las  Cruces,  N.  Mex.  223-328. 

16.  Amlicke,  B.B.,  and  Coleman,  I.W.  (1984)  High  Resolution  Wind  (IIRWJ 

Model,  MRC  /  W  DC -R -089,  Mission  Research  Corporation,  Alexandria,  Va. 

17.  I.anczos,  C.  (1970)  The  Variational  Principles  of  Mechanics  (4th  ed),  L.  of 

Toronto  Press,  Toronto,  pp.  106-110. 


Table  1.  Description  of  Windflow  Models  Considered 


Model 

Descript  ion 

Disadvantages 

U.S.  Army  Atmospheric 
Sciences  Lab  (ASL)^ 

A  two-dimensional  surface 

layer  model  suitable  to  small 
grid  spacings  and  complex 

terrain.  Considers  anabatic 

and  katabatic -like  flows  influ¬ 
enced  by  static  stability. 

Model  is  run  on  a  UNIVAC  1108 

computer. 

Considers  surface 

layer  only.  Extra¬ 
polation  into  mixed 

layer  is  not  recom¬ 
mended. 

Oak  Ridge,  Tenn.  ^ 

A  two-dimensional  model 

which  employs  variational 
analysis  technique  for  use  in 
conjunction  with  a  movable 
fine  mesh,  allowing  a 
Lagrangian  solution  of  the 
advection-diffusion  equation 

(ADE)  to  be  calculated. 

If  shear  or  turbu¬ 
lence  on  a  scale 

larger  than  the  scale 

of  the  mesh  but 

smaller  than  the 

scale  of  the  full 

mesh  exists,  the 

mesh  may  become 

distorted  in  time, 

leading  to  numerical 

errors  in  the  solution. 

MASCON  (Lawrence 
Livermore  Lab^) 

A  two-dimensional  mass- 

conservation  model  which 

allows  the  mixed-layer 
depth  to  vary  in  time  and 

space.  Horizontal  flux 
fields  are  realistically 
produced  through  this 
method.  Channeling  ef¬ 
fects  of  topography  are  also 
simulated  through  use  of 
terrain  (no-flow  through 
boundaries)  and  temporally 
varying  inversion  heights. 

Horizontal  resolu¬ 
tions  below  500  m 

are  not  advised  be¬ 
cause  of  problems 
associated  with  spe¬ 
cification  of  terrain 

(obstacle  cells  re¬ 
present  terrain). 

The  model  currently 

runs  on  a  CDC  7600 

and  takes  between  50 

and  70  iterations  to 

converge  to  a 

steady -state -solu¬ 
tion.  This  means 

model  will  take  a 

long  time  to  run  on 
a  minicomputer  in 

present  form.  Wind 
velocity  problems 

exist  for  regions 
where  the  depth  of 
the  mixed  layer  is 

shallow  because  of 

high  terrain  and  low 
inversion  height. 

Table  1  (cont.  ).  Description  of  Windflow  Models  Considered 


Model 

Description 

Disadvantages 

University  of 
Washington"7 

A  two-dimensional  primitive 
equation  (PE)  model  which 
can  use  operational  tempera¬ 
ture  and  wind  soundings  for 
input  data  (usually  at  850  mb). 
Model  is  diabatic  and  adapt¬ 
able  to  run  on  PRIME  400. 

Cyber  750.  and  CRAY  1 

computers. 

Horizontal  grid  re¬ 
solutions  below 

3km  are  not  recom¬ 
mended  because  of 

the  breakdown  of  the 

hydrostatic  approxi¬ 
mation  below  dx  «■<- 
3  km  and  Az  ^  1  km®'  ® 

Lavoie10'11 

A  two-dimensional  mixed- 

layer  model  which  employs 
variable  horizontal  grid 
spacing  and  parameterization 
of  clouds  and  precipitation. 

Model  takes  280s  CPU  time 

on  IBM  360-67. 

Model  assumes  a 

well-mixed  planet¬ 
ary  boundary  layer 
(PEL).  Hydrostatic 
assumption  breaks 

down  for  horizontal 

resolutions  needed 

in  our  case. 

WINDS  (LSD A  and 

Colorado  State 

University^) 

An  inexpensive  two-dimen¬ 
sional  model  which  produces 
reasonably  accurate  wind 
speeds  and  directions  while 
incorporating  effects  of  sur¬ 
face  heating  and  topography. 

A  mass-conservation  model, 

it  uses  surface  temperatures, 
pressures,  and  roughness 
lengths  to  modify  potential 
flow  and  employs  additions  of 
divergence  and  vorticity  fields 

to  the  basic  flow  as  well. 

Output  only  yields 

wind  flow  between  1 

and  2  km  above 

ground  level  (ACL), 
so  it  is  impractical 

for  use  as  a  diagnos¬ 
tic  model  for  surface 

wind  flow. 
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where  is  the  total  domain  residual  as  computed  from  Eq.  (1),  and  R^  is  the  in¬ 
dividual  residual  at  each  grid  point  i,  j.  A  number  of  relaxation  sweeps  are  then 
performed  in  the  model  to  compute  a  velocity  correction  proportional  to  the  steep¬ 
est  descent  vector  in  the  multi-dimensional  velocity  space,  written  in  Eq.  (3): 


n*  *  -  (dR.Jdvhir  £  <3R tJdvh2]112 

H  1J  Li,j,k  x3  XJ 


where  k  -  1,  2  and  denotes  the  u,  v  wind  components,  n  is  the  nondimensional  cor¬ 
rection  factor,  and  v  is  the  velocity.  The  velocity  correction  is  given  by  the  fol¬ 
lowing  relation: 


dv^  =  0.03  vq  /n  x  M  n^. 


where  4v.^  is  the  velocity  correction  at  i,  j,  vq  is  the  initially  uniform  velocity, 
and  N,  M  are  the  grid  dimensions.  The  relaxation  sweeps  and  windfield  adjust¬ 
ments  are  performed  until  the  residual  of  the  integral  in  Eq.  (1)  begins  to  "con¬ 
verge"  to  a  minimum  value. 

There  arc  several  reasons  why  a  variational  approach  is  appropriate  for  this 

4 

modeling  application.  Ball  and  Johnson  stated  that  the  advection  terms  of  the 
momentum  equation  are  dominant  at  this  scale,  and  it  is  obvious  that  thermal  stra¬ 
tification  is  important  in  the  forcing  of  topographic  flows.  Variational  techniques 

have  been  used  to  derive  the  equation  of  motion  for  inviscid,  incompressible 
4  18 

flow.  ’  Considering  the  points  mentioned  above  and  the  small  space  scale  in¬ 
volved,  it  makes  sense  to  use  a  variational  approach  to  the  simulation  of  the  sur¬ 
face  wind  field  over  complex  terrain  considering  advection,  buoyancy*,  detailed 
topography,  and  mass  conservation  (a  consequence  of  the  incompressibility  as¬ 
sumption). 

The  model  domain  covers  a  limited  area  (5  by  5  km)  located  in  the  Kisatchie 
National  Forest  on  the  Fort  Polk  Military  Reservation,  about  25  km  ESE  of 


♦The  "buoyancy"  referred  to  in  this  context  [expressed  in  Eq.  (7)]  is  not  a 
strictly  vertical  variation  because  of  the  terrain -following  coordinate  system.  A 
horizontal  temperature  gradient  is  introduced  in  this  system;  thus,  a  horizontal 
pressure  gradient  force  appears  which  is  the  actual  driving  mechanism  involved 
in  slope  flows. 

18.  Dutton,  J.A.  (1976)  The  Ceaseless  Wind.  McGraw-Hill,  New  York,  pp.  435- 
439. 
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Leesville,  La.  The  area  consists  of  rolling  hills  with  tall  vegetation.  ^  The 
model  makes  use  of  limited  meteorological  input  data  (one  observation  point,  using 
the  surface  wind,  and  potential  temperatures  and  geopotential  heights  at  two  levels 
above  the  surface).  The  model  employs  detailed  terrain  and  vegetation  data  (over 
a  51-by-51  grid  with  a  grid  spacing  of  100  m),  the  latter  used  to  compute  surface 
roughnesses  and  modify  the  terrain  elevations.  The  model  uses  a  staggered  grid 
to  solve  for  wind  velocity  at  grid  points  between  the  terrain  points  so  that  terrain 
slopes  can  be  computed. 

These  model  characteristics  are  essentially  the  same  in  both  ASL  and  AFGL 
versions.  The  modifications  made  to  the  ASL  model  are  described  in  the  next 
section. 

2.2  Modifications  Made  to  the  Model  at  AFGL 

The  ASL  windflow  model  had  several  computational  algorithms  that  we  exam¬ 
ined  in  great  detail.  The  first  of  these  was  the  determination  of  surface  pressure 
py  and  "reference"  potential  temperature  0^,  the  latter  used  as  a  comparison 
against  a  calculated  surface  potential  temperature  for  buoyancy  determination. 

In  the  ASL  version  of  the  model,  both  pg  and  6q  are  calculated  through  a  linear 
extrapolation  technique  that  uses  geopotential  heights  and  potential  temperatures 
at  two  levels  above  the  surface  (850  and  700  mb). 

We  addressed  the  pg  calculation  first.  We  noted  that  using  a  linear  extrapola¬ 
tion  technique  for  a  nonlinearly  varying  quantity  such  as  p  can  lead  to  errors  in 
the  determination  of  pg  and  to  subsequent  errors  in  the  calculations.  We  deter¬ 
mined  that  if  a  linear  extrapolation  is  used,  it  should  encompass  as  small  a  Ap  as 
possible  and  be  done  close  to  the  surface.  We  chose  to  do  the  pg  calculations  by 
using  geopotential  heights  at  950  and  900  mb,  since  these  levels  are  close  to  the 
surface  over  Fort  Polk  and  yield  more  realistic  values  of  py,  as  illustrated  in  the 
example  shown  in  Table  2. 

When  we  addressed  the  linear  extrapolation  of  0Q,  it  became  clear  that  the 
reference  6()  calculated  this  way  does  not  represent  a  true  dry  adiabatic  lapse 
rate  from  850  mb  down  to  the  surface.  Thus,  any  buoyancy  calculation  using  this 
0o  value  is  not  an  accurate  estimation  of  the  actual  stability  conditions.  Moreover, 
similar  errors  to  those  in  pg  result  from  extrapolation  using  0  values  at  850  and 
700  mt>.  Since  the  buoyancy  determination  in  the  model  determines  the  stability 
regime  for  the  surface  layer,  we  felt  that  the  950-mb  potential  temperature  would 
tie  more  representative  of  the  average  0  in  the  mixed  layer,  and  that  comparing 

19.  Sadeh,  W.X.,  Law,  F.  W.  .  Marlatt,  W.  K. ,  Fox,  I).  CL.  Dietrich.  D.  L. ,  and 
C'ionco,  R.M.  (1982)  A  Survey  of  Micromoteorological  Parameters  Within 
a  Forest  C  anopy  at  Fort  Polk,  Louisiana,  AsT.-d? -0100 -I,  .  Army 

Atmospheric  Sciences  Laboratory,  White  Sands  Missile  Range,  N.  Mex. 
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Table  2.  Sample  Calculations  of  ps  Using  Linear  Extrapolation  From  950,  900  mb 
as  Opposed  to  850,  700  mb 


Pu  '  PL 

Standard  Formula:  p  *  c - r —  (h  -  hT  )  +  p7  , 

s  n,,  -  hr  s  L  L 


where  s  =  surface,  u  =  upper  level,  and  L  =  lower  level. 

Let  h  =  75  m  (average  elevation  over  Fort  Polk),  and  use  climatological  h 


values  for  January*: 


h9S0  =  599  m 

C-“«- 

Using  950,  900  mb: 
Using  850,  700  mb 


h850  =  1520  m 

h700  =  3115  m 

p  =  1008.  5  mb 
*s 

p  =  985.  9  mb 
rs 


|  Using  hyposometric  eq  (T  =  283.  15  K):  pg  =  1012. 0  mb 

♦Climatological  values  for  h,  T  are  taken  from  upper-air  data  for  Lake  Charles, 


this  0Q  with  Q  would  give  a  more  meaningful  estimation  of  the  surface-layer  sta¬ 
bility  than  comparison  to  a  linearly  extrapolated  using  0  from  850  and  700  mb. 
This  point  is  illustrated  in  Figure  1.  Once  this  change  in  the  $  determination 
has  been  made,  it  is  no  longer  necessary  to  adjust  the  buoyancy  in  the  model  on 
the  basis  of  a  dQ/dZ  computation  (as  is  done  in  the  ASL  version:  see  Table  2). 

Another  change  made  to  the  ASL  version  of  the  model  was  the  input  for  deter¬ 
mining  $s  and  buoyancy.  In  the  ASL  version,  a  surface  heating  rate  A Q  (based  on 
first  law  of  thermodynamics  considerations)  is  input  into  the  model,  and,  from  it, 
values  of  0g,  T  (surface  air  temperature),  and  buoyancy  are  determined.  We 
believed  that,  in  an  operational  scenario,  it  would  be  difficult  to  estimate  AQ,  but 
Tg  is  readily  available.  In  the  AFGL  version,  Tg  is  an  input  parameter,  and  it  is 
used  to  calculate  a  field  based  on  the  pg  value  at  each  grid  point.  This 
field,  when  compared  with  the  0^  field  (the  latter  being  a  constant  field  equal  to 
^950*'  £*ves  information  about  the  buoyancy  as  a  function  of  terrain  elevation.  It 
should  be  noted  here  that  the  assumption  of  a  constant  $  field  under  daytime* 
conditions  is  reasonable  over  a  region  such  as  Fort  Polk,  whose  terrain  elevations 
only  vary  between  60  and  95  m  above  sea  level.  However,  over  a  more  mountain¬ 
ous  region,  this  assumption  should  be  revised. 

*At  night,  choice  of  an  appropriate  can  be  difficult  because  of  the  formation  of 
the  nocturnal  inversion.  This  problem  is  then  complicated  by  the  presence  of 
mountains. 


:ltV 


q>  (A  S  L  version) 


a»o=q3  (A.F6  L.  version) 


BUOYANCY  (A  FG  L  version) 


BUOYANCY  (A.SL  version) 

Figure  1.  Schematic  Diagram  Illustrating  Differences  in  Buoyancy  Calculation  for 
the  ASL  and  AFGL  Versions  of  the  Model.  The  thick  solid  line  denotes  a  potential 
temperature  sounding  for  the  unstable  case.  The  dashed  line  is  the  calculation  of 
0o  by  linear  extrapolation  using  0  at  850  and  700  mb.  HPBL  refers  to  the  plane¬ 
tary  boundary  layer  depth,  defined  here  as  the  height  at  which  a  stable  lapse  rate 
is  first  encountered.  In  this  example,  =  8950  (representative  of  the  mean  po¬ 
tential  temperature  in  the  mixed  layer)  gives  a  much  more  realistic  measure  of 
the  buoyancy  associated  with  unstable  conditions  than  the  linearly  extrapolated  0O 
does 


Finally,  a  word  should  be  said  about  calculation  of  the  bulk  Richardson  num¬ 
ber  Ri  in  the  model.  The  ASL  version  of  the  model  includes  a  "curvature"  Rich¬ 
ardson  number  based  on  the  effects  of  flow  curvature  on  stability  as  documented 
in  a  paper  by  Bradshaw.  However,  Bradshaw's  work  was  based  on  physical 
modeling  studies  of  boundary  layers  on  curved  walls,  and  implies  a  relationship 
between  temperature  fluctuations  (producing  buoyancy  fluctuations)  and  u-wind 
component  fluctuations  (producing  Coriolis  force  fluctuations).  A  scale  analysis 

2  j 

of  mesoscale  atmospheric  motions  performed  by  Atkinson  shows  that,  for  time 

and  space  scales  of  the  order  that  we  are  working  with  for  the  Fort  Polk  experi- 
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ments,  the  Coriolis  force  is  relatively  unimportant.  Also,  Walmsley  says  that 
in  the  surface  layer,  the  eddy  Ekman  number  is  usually  much  greater  than  1,  al¬ 
lowing  the  Coriolis  force  to  be  neglected.  Thus,  we  felt  that  a  curvature  Ri  as  in¬ 
cluded  in  the  ASL  model  was  inappropriate  for  our  work.  In  several  preliminary 


20.  Bradshaw,  P.  (1969)  The  analogy  between  streamline  curvature  and  buoyancy 

in  turbulent  shear  flow,  J.  Fluid  Mech,  36.,  Part  1:177-191. 

21.  Atkinson,  B.  W.  (1981)  Mesoscale  Atmospheric  Circulations.  Academic 

Press,  London,  pp.  13-20. 


Table  3.  Comparison  of  ASL  and  AFGL  Versions  of  Model 


ASL  Version 

AFGL  Version 

Uses  850,  700  mb  h,  6  for 

Uses  950,  900  mb  h  for  linear 

linear  extrapolation  of  pg,  0Q • 

extrapolation  of  pg. 

Uses  linearly  extrapolated  0Q 

Takes  #ggQ  and  uses  it  as 

as  a  reference  potential  tem- 

reference  mixed-layer  potential 

perature  against  which  0g  is 

temperature  0Q  against  which 

compared  to  determine  value 

#g  is  compared  to  determine 

of  buoyancy. 

value  of  buoyancy. 

Uses  curvature  Ri  to  modify 

Bulk  Ri  alone  is  used  to  deter- 

bulk  Ri  to  allow  for  effects  of 

flow  curvature  on  stability. 

mine  stability. 

Uses  surface  heating  4Q  as 

Uses  observed  Tg  as  input  and 

input  parameter  for  calcula- 

then  calculates  0  from  T 
s  s 

tion  of  6  and  T  . 

and  p  . 

s  s 

Does  an  adjustment  to  0g, 

Since  buoyancy  is  calculated 

buoyancy,  and  Ri to  account 

based  on  comparison  with  poten- 

for  lapse  rate  d6ld z  using 

tial  temperature  at  950  mb,  no 

power  law  type  formulation. 

adjustment  of  0g,  buoyancy,  or 
Ri  is  necessary. 

experiments  with  and  without  this  eurvature  Hi,  our  hypothesis  was  confirmed  as 
the  flow  fields  were  essentially  the  same  in  both  experiments  (figures  not  shown). 

Tlie  major  differences  between  the  ASI.  and  AFGL  versions  of  the  windflow 
model  are  summarized  in  Table  3, 


3.  MODEL  EXPERIMENTS 


This  section  describes  the  results  of  our  model  evaluation  phase.  Analysis 
includes  detailed  examinations  of  kinetic  energy  in  the  model,  and  the  power  law 
assumption  as  applied  to  the  wind  profiles  in  the  surface  layer.  We  then  present 
the  results  of  our  sensitivity  experiments. 


For  the  model  experiments,  we  used  simulated  data  based  on  surface  clima¬ 
tology  for  Fort  Polk  and  upper  air  climatology  for  nearby  Lake  Charles.  The 
simulated  data  were  used  to  construct  soundings  representative  of  the  atmosphere 
for  typical  January  mornings  (stable  conditions)  and  July  afternoons  (unstable  con¬ 
ditions).  The  July  afternoon  soundings  were  also  used  to  construct  experiments 
simulating  neutral  conditions. 

3.1  Fort  Polk  Topography 

19 

A  detailed  description  of  the  Fort  Polk  topography  appears  in  Sadeh  et  al. 

A  chart  showing  the  terrain  and  vegetation  over  the  model  domain  is  shown  in  Fig¬ 
ure  2.  Notice  from  Figure  2  that  there  are  primarily  two  lowland  areas,  one  run¬ 
ning  northwest  to  southeast  across  the  center  of  the  domain  (corresponding  to  Six- 
mile  Creek),  and  another  in  the  southwest  corner  of  the  domain.  A  ridge  runs 
parallel  to  the  Sixmile  Creek  lowland  area,  and  high  terrain  also  appears  at  the 

northern  end.  The  tallest  vegetation  over  the  domain  (coniferous  trees)  averages 
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between  21.3  and  22.5  m  high  with  an  average  density  of  1  tree/4.  22  m  . 

Detailed  terrain  and  vegetation  data  were  supplied  to  us  for  this  area. 

Although  the  terrain  elevations  over  this  region  only  varied  by  10s  of  m,  the  pre¬ 
sence  of  tall  vegetation  could  enhance  the  complexity  of  the  topography  consider¬ 
ably.  For  instance,  a  terrain  slope  from  lowland  marshes  to  pine  tree-capped 
hills  could  be  as  great  as  40  to  50  m  over  a  distance  of  only  several  hundred  me¬ 
ters.  Additionally,  the  presence  of  nonhomogeneous  vegetation  could  produce  lo¬ 
calized  "low  spots,  "  which  could  focus  the  locations  of  cold  pools  of  air  during 
nocturnal  drainage  episodes.  The  model  considers  the  effects  of  vegetation  in  two 
ways:  (1)  using  vegetation  heights  to  modify  the  terrain  elevation,  and  (2)  through 
the  surface  roughness.  Lq.  (5)  illustrates  the  terrain  modification  by  vegetation: 

h'  h  •  0.  7  h  ,  (5) 

o  veg 

where  h'  is  the  modified  terrain  height,  is  the  original  terrain  height,  and  h 

is  the  vegetation  height.  Lq.  (5)  describes  the  concept  of  "displacement  length,  " 

which  is  given  as  being  between  70  and  80  percent  of  hveg  and  is  used  to  define  a 

height  h'  to  which  the  logarithmic  wind  profile  can  be  extrapolated  downward  from 
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above  (see  Figure  6.2  in  Panofsky  and  Dutton  ).  The  0.7  hveg  term  in  Kq.  (5) 
agrees  well  with  the  0.  68  hveg  form  proposed  by  Monteith.  Figure  3  illustrates 

22.  Panofsky,  H.A.,  and  Dutton,  J.A,  (1984)  Atmospheric  Turbulence:  Models 
and  Methods  for  Kngineering  Applications,  Wiley,  New  York,  pp.  124-125. 

2.3.  Monteith,  J.  L.  (1973)  Principles  of  Lnvironmental  Physics.  Edward  Arnold 
Ltd.,  London,  pp.  88-90. 


the  topography  over  the  domain  when  Eq.  (5)  is  applied  to  the  terrain  elevations. 

An  examination  of  Figure  3  reveals  the  appearance  of  several  "low  spots"  over  the 
domain,  with  elevations  h'  less  than  80  m.  Also,  several  "peaks"  have  appeared 
corresponding  to  tall  pines  over  high  ground,  with  elevations  over  110  m.  Addi¬ 
tionally,  the  nonhomogeneity  of  the  terrain  has  been  increased  through  the  addition 
of  vegetation.  A  comparison  of  Figures  2  and  3  reveals  that  t tie  addition  of  vege¬ 
tation  information  to  the  terrain  heights  alters  the  character  of  the  topography 
over  the  Fort  Polk  region. 

The  second  way  of  considering  the  effects  of  vegetation  on  the  topography  is 

through  the  surface  roughness.  Many  formulations  exist  for  considering  roughness 
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as  a  function  of  vegetation  height  (see  Monteith,  '  Sellers,  and  Lettau  for  dis¬ 
cussion  of  these  formulations),  in  the  model,  roughness  Zq  is  computed  for  each 
grid  point  by  the  following  formula: 

X  =  0.  14  (h  *  0.  1).  (6) 

o  veg  '  ' 

Table  4  shows  the  Z()  values  computed  from  Kq.  (6)  for  various  vegetation  types 
found  over  the  domain.  The  values  of  '/.  in  Table  4  compare  favorably  with  those 


Table  4.  Surface  Roughness  Values  Derived  From  Vegetation 
Heights 


Vegetation  Type 

h  (m) 

veg 

Zo  (cm) 

Deciduous  and 

Coniferous  Trees 

21.  3 

299.  7 

Coniferous  Trees 

20.  5 

288.  4 

Small  Coniferous  Trees 

14.  5 

204.  4 

Marshland 

0.  9 

14.  1 

<  learing 

0.01 

1.5 

24.  Sellers,  VV.l).  (1965)  Physical  Climatology,  1.  of  Chicago  Press,  Chicago, 

pp  150-151. 

25.  I.ettau,  11.  (1969)  Note  on  aerodynamic  roughness-parameter  estimation  on 

the  bases  of  roughness-element  description,  J.  Appl.  Meteorol.  8:828-832. 
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obtained  using  Tanner  and  Pelton's  logarithmic  formula,  setting  a  -  -0.883  and 
b  =  0.  997. 

3.2  Analysis  of  Model  Parameterization* 

3.  2.  1  WIND  PROFILE  SENSITIVITIES  TO  STABILITY  AND  ROUGHNESS 

The  windflow  model  determines  stability  through  a  buoyancy  calculation  given 
by  the  following  equation: 


B  =  g 


(7) 


where  g  is  the  gravitational  acceleration,  -  0g5Q  (see  Section  2.2),  and  &s  is 
the  surface  potential  temperature.  This  value  of  buoyancy  is  then  inserted  into 
the  bulk  Richardson  number  calculation,  given  by  Eq.  (8)  below: 


Ri  =  B  Z/u2. 


(8) 


where  Z  =  ^GRID  +  Z0  ^GRID  a  reference  height  -  10  m,  in  our  case),  and  u  is 
the  windspeed  at  Z.  For  windspeeds  u  ^  0.  1  ms'1,  the  bulk  Ri  formula  is 

Ri  =  100  BZ.  (9) 


This  form  of  Ri  produces  large  magnitudes  of  Ri,  which  explains  in  part  the  reason 
that  the  model  can  reproduce  topographically  driven  flows  given  initially  calm 
wind  conditions  (see  Section  3.3).  Determination  of  Ri  defines  the  stability  class¬ 
ification  of  the  surface  layer.  Seven  stability  "classes"  are  used  in  the  model, 
based  on  Ri  values  from  -5,0  through  0.25.  These  stability  classes  are  used  to 
calculate  values  of  the  Monin-Obukhov  similarity  parameters,  which  help  to  deter¬ 
mine  the  value  of  the  wind  profile  exponent  n,  given  by  the  following  relations: 

{ <t>ml  ?  n  Z./Z0  -  lU  ']  ,  Ri  <  0.0 

n  -  ■  n  Z/Z0  +  4.  7  (x  -  xQ)]  ,  0.  0  Ri  <  0.  25  (10) 

{  1  ,  Ri  ^  0.  25. 

where  0  is  the  nondimensional  wind  shear,  l^'is  a  Monin-Obukhov  function,  and 


26.  Tanner,  C.  B.  ,  and  Belton,  W.  L.  (1960)  Potential  evapotranspiration  esti¬ 
mates  by  the  approximate  energy  balance  method  of  Penman,  J.  Geophys. 
Res.  65:3391-3413.  - - 


x.  xo  are  proportional  to  (Z/L)  and  (Zq/l)  11  \  respectively  [see  Eqs.  (2.30 
a-d)]  in  Amlicke  and  Coleman. 

The  model  makes  use  of  the  Husinger  et  al27  wind  profiles  by  assuming  that 
they  can  be  applied  locally  at  each  grid  point.4  The  wind  profile  exponent  is  used 
in  a  power  law  formation  to  evaluate  the  wind  at  the  computation  height  Z  (top  of 
the  flux  box).  The  values  of  used  in  the  model  are  the  same  as  those  of 
Husinger  et  al  and  are  not  directly  influenced  by  different  Zo  values  (see  Figure  4) 


<£mvs  Ri  •  o'l  Z0 


Figure  4.  Graph  .>f  <t>m  (Non-Dimensional  Wind  Shear)  as  a  Function  of  the  Hulk 
Richardson  Number  Ri.  for  All  Values  of  Surface  Roughness  ZQ 


27.  Husinger,  J.A.,  Wvngaard,  .T.C.,  bumi,  Y .  .  and  Hradlev,  F.F.  (1971) 
Flux-profile  relationships  in  the  atmospheric  boundary  layer.  J.  Atmos 
Sci.  28:181-189.  '  - — 1 - - 


the  value  of  n  is  sensitive  to  the  roughness,  as  illustrated  by  the  graph 
5.  Figure  5  shows  the  value  of  n  increases  with  increasing  Z.^.  Thus, 
ble  for  roughness  effects  to  feed  back  to  the  windfield  during  the  model 
l  and  alter  the  stability  of  the  surface  layer. 


3.2.  2  ANALYSIS  OF  MOIJKl.  KINFTIC  FNEROY 


The  basic  model  Kq.  (1)  describes  the  adjustment  of  an  initially  uniform  wind- 
field  (based  on  a  single  input  observation)  to  conform  with  buoyancy  and  topography. 
During  this  adjustment,  the  integral  residual  R  must  be  continually  decreasing 
until  a  quasi-steady  state  balance  is  achieved  and  the  value  of  R  reaches  a  mini¬ 
mum.  Since  the  buoyancy  essentially  remains  the  same  throughout  the  model  re¬ 
laxation,  it  follows  that  a  considerable  adjustment  in  the-  total  domain  kinetic  en¬ 
ergy  KF  takes  place  during  this  time. 

We  studied  the  changes  in  total  KK  during  the  mod  1  integrations  in  detail, 
examining  sensitivity  to  stability  (that  is,  buoyancy  I,  initial  wind  -peed,  and  num¬ 
ber  of  model  iterations.  Our  results  can  be  summarized  bv  the  following  state¬ 
ments: 

(1)  KK  is  sensitive  to  the  magnitude  of  the  buoyancy  in  the  model.  Final  KK 
is  higher  for  larger  magnitudes  of  buoyancy.  KK  decreases  through  the  integra¬ 
tion  for  the  neutral  case. 

(2)  KK  is  sensitive  to  the  initial  wind  speed.  When  initially  low  (2  ms  *)  wind 
speeds  are  used,  the  KK  tends  to  increase  through  the  integration.  When  the  ini¬ 
tial  wind  speeds  are  high  (10  ms  *),  KK  decreases  through  the  integration. 

(3)  The  model  wind  field  essentially  reaches  quasi-steady  state  after  60  iter¬ 
ations. 

We  will  now  address  each  of  these  results  individually. 

The  first  result  makes  sense  when  one  examines  Kq.  (1),  since  the  winds  in 
the  acceleration  term  would  have  to  tie  stronger  to  adjust  to  a  stronger  value  of  B 
in  the  integral.  This  point  is  further  addressed  in  Section  3.3.  The  sensitivity  of 
the  model  KK  to  buoyancy  is  illustrated  by  the  graphs  in  Figure  6,  which  show  the 
change  in  KK  as  a  function  of  the  number  of  iterations  NR.  Note  that  in  the  exper¬ 
iment  with  B  =  -0.06  (slightly  unstable),  the  KK  actually  decreases  slightly  during 
the  model  integration.  This  occurs  because  of  the  small  magnitude  of  buoyancy 
used  in  that  experiment.  The  experiment  with  B  =  -0.  18  (moderately  unstable), 
which  also  shows  a  decrease  in  KK,  is  one  in  which  an  initial  wind  speed  of 
10  ms  1  was  used.  An  initial  windspeed  of  10  ms  *  was  also  used  for  the  neutral 
case  (B  0  curve  in  Figure  6).  The  results  from  Figure  6  suggest  that  there 
exists  some  intermediate  wind  speed  between  2  and  10  ms  *  for  which  KK  would 
remain  fairly  constant  through  the  model  integration.  As  far  as  the  third  result 
is  concerned,  notice  that  Figure  6  shows  the  KK  changes  for  four  experiments 
which  were  integrated  to  NR  200.  A  total  of  eight  experiments  were  actually 
run  with  NR  =  200,  and  the  results  of  these  experiments  indicated  that,  at  NR  -  60, 
KE  reached  roughly  65  percent  of  its  value  at  NR  =  200.  This  percentage  also  ap¬ 
pears  to  he  sensitive  to  the  magnitude  of  the  buoyancy,  as  shown  in  Figure  6, 


However,  a  comparison  of  the  wind  fields  produced  at  NR  =  60  with  those  at  NR  » 
200  (figures  not  shown)  reveals  that  the  flow  patterns  fully  develop  in  the  model  by 
NR  =  60. 

3.3  WindfloH  fharai'lrrislirii  I'nder  Various  Stability  Conditions 

This  section  documents  our  results  from  the  model  sensitivity  experiments 
with  and  without  vegetation.  The  formulation  of  soundings  simulating  stable,  un¬ 
stable,  and  neutral  conditions  was  described  at  the  beginning  of  Section  3.  The 
climatological  soundings  for  the  two  stable  cases  are  shown  in  Figure  7a.  Two 
soundings  each  were  formulated  for  this  and  the  unstable  case  to  test  the  model 
sensitivity  to  the  magnitude  of  the  buoyancy.  Similar  soundings  for  the  unstable 
cases  are  shown  in  Figure  7b.  The  neutral  case  sounding  is  shown  in  Figure  7c. 

Climatology  for  Fort  Polk  indicated  that  the  preferred  wind  direction  and 
speed  during  January  mornings  (0600-0800  LST)  was  from  the  north  at  2.  3  ms  * 
(although  calm  actually  has  the  highest  frequency  for  this  time).  The  southerly 
direction  was  the  second  most  frequent  direction,  and  this  direction  was  chosen 
for  the  sensitivity  tests  to  allow  direct  comparisons  with  the  unstable  cases.  For 
July  afternoons  (1500-1700  LST),  south  was  the  preferred  direction  with  a  mean 
speed  of  3.0  ms  *.  Calm  had  the  second  highest  frequency  at  this  time. 

Table  5  gives  the  input  information  used  to  initialize  the  model  for  the  eight 
experiments  that  tested  model  sensitivity  to  stability  and  wind  conditions.  It 
should  be  noted  here  that  experiments  were  run  for  the  model  using  other  wind  di¬ 
rections,  and  the  results  were  consistent  with  the  model  physics  and  results  using 
a  southerly  wind. 

The  model  domain  for  the  experiments  with  vegetation  appeared  in  Figure  3, 
and  the  domain  for  the  experiments  without  vegetation  appeared  in  Figure  2.  The 
model  runs  made  using  the  vegetation  information  employed  horizontally  varying 
roughness.  The  experiments  with  no  vegetation  used  a  horizontally  uniform  rough¬ 
ness  of  14  cm,  corresponding  to  h  =  0.  9  m  (rolling  hills  with  tall  grasses). 

veg 

This  value  of  Zq  was  chosen  in  order  not  to  be  on  the  lowest  wind  exponent  curve 
(see  Figure  5),  yet  realistically  simulate  a  lack  of  tall  vegetation  in  the  area. 

A  problem  exists  when  trying  to  interpret  results  from  experiments  using  no 
vegetation  and  uniform  Zq.  The  problem  is  how  to  isolate  the  effects  on  the  flow 
field  of  vegetation  alone  from  those  caused  by  variable  roughness  alone.  To  solve 
this  problem,  model  runs  were  made  where  the  effects  of  vegetation  on  the  terrain 
heights  were  maintained,  but  a  uniform  of  14  cm  was  used.  These  results 
could  he  compared  with  results  from  model  runs  using  vegetation  and  variable  Z 
to  isolate  the  effects  on  the  windfield  of  variable  roughness  alone.  These  compar¬ 
isons  revealed  that  the  effects  of  horizontally  varying  roughness  are  minimal. 
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Experiment  No.  ♦* 

Wind  Dir 

Speed  (ms'*) 

T  g  (°  C  ) 

1  (stable) 

S 

2 

7.2 

2  (very  stable) 

S 

2 

-0.4 

3  (slightly  unstable) 

s 

2 

30.  6 

4  (very  unstable) 

s 

2 

34.4 

5  (very  unstable) 

s 

10 

34.4 

6  (neutral) 

s 

10 

28.  7 

7  (very  stable) 

Calm 

-0.4 

8  (very  unstable) 

Calm 

34.4 

B*(ms'^) 


♦Values  of  B  based  on  6g  values  calculated  using  pg  (Jan)  =  1007.  4  mb  and  pg  (July) 
*  1004.0  mb. 

♦♦Experiments  using  no  vegetation,  with  uniform  roughness  of  14  cm  are  numbered 
the  same  as  below  with  a  suffix  "a.  " 
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The  flow  patterns  remained  essentially  the  same,  with  the  upslope  and  drainage 
flows  nearly  identical  in  experiments  with  and  without  nonuniform  roughness  (fig¬ 
ures  not  shown).  One  curious  result  of  the  comparison  was  that  the  total  KE  over 
the  domain  averaged  about  13  percent  less  for  the  uniform  Zq  cases  than  for  the 
varying  Zq  cases.  We  believe  these  differences  to  be  rather  small,  considering 
the  total  grid  size  involved  (50  by  50  grid;  2500  grid  points).  These  results  sug¬ 
gest  that  the  model  runs  using  vegetation  and  nonuniform  roughness  can  be  di¬ 
rectly  compared  with  experiments  using  no  vegetation  and  uniform  roughness,  be¬ 
cause  the  effects  of  variable  roughness  are  second-order. 

3.3.1  STABLE  CONDITIONS  WITH  VEGETATION 

The  stable  surface  layer  was  simulated  in  experiments  1,  2,  and  7.  The  input 
data  for  these  experiments  are  shown  in  Table  5.  In  experiments  1  and  2,  the  in¬ 
put  wind  speed  of  2  ms  *  is  close  to  the  climatological  value  of  3.5  ms'1  for  this’ 
direction  in  January. 

Figure  8a  shows  the  windfield  for  experiment  1.  Notice  that  there  are  several 
convergence  zones  corresponding  to  cold  air  drainage  into  low  spots  on  the  domain. 
The  model  is  capable  of  simulating  cold  air  drainage  off  the  higher  terrain  and 
flow  diversion  around  terrain  obstacles  such  as  hills  and  ridges.  When  the  buoy¬ 
ancy  is  increased  (experiment  2;  Figure  8b),  the  windfield  is  essentially  the  same 
except  for  some  weakening  of  the  flow  over  the  level  terrain  areas  and  some  en¬ 
hancement  of  the  slope  winds  in  the  drainage  areas.  The  enhancement  of  the  slope 
flow  occurs  because  of  the  model  response  to  buoyancy  and  topography.  Recall 
from  Section  2  that  the  surface  layer  is  divided  into  flux  boxes  of  thickness  Z  and 
that  these  boxes  are  distorted  to  conform  with  the  terrain  slopes.  Also  recall 
from  Section  3.  2.  2  that  model  KE  was  larger  for  experiments  with  larger  magni¬ 
tudes  of  buoyancy.  At  this  point  in  our  evaluation,  it  is  too  early  to  verify  this 

slope-flow  enhancement  against  observations,  although  these  results  are  inconsis- 

28 

tent  with  the  numerical  results  of  McNider  and  Pielke  and  analytical  results  of 

2Q 

McNider. 

The  model  is  also  capable  of  generating  cold  air  drainage  flows  given  initially 
calm  wind  conditions.  This  result  is  illustrated  in  Figure  8c,  which  shows  the 
model  windfield  for  experiment  7.  Notice  that  the  drainage  is  focused  into  low- 
lying  areas  (lowland  marshes  along  Sixmile  Creek),  and  that  the  strongest  drain¬ 
age  occurs  over  the  steepest  terrain  slopes.  Again,  these  results  are  consistent 
with  the  constraint  Eq,  (1). 


28.  McNider,  R.T.,  and  Pielke,  R.A,  (1984)  Numerical  simulation  of  slope  and 

mountain  flows,  J.  Climate  and  Appl,  Meteorol.  23:1441-1453. 

29.  McNider,  R.T.  (1982)  A  note  on  velocity  fluctuations  in  drainage  flows, 

J.  Atmos.  Sci,  39:1658-1660. 
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Figure  8a.  Model  Windfield  for  Fxperiment  1  Showing  Wind  Velocity  Vectors  for 
Stable  C’a.se  With  Vegetation  (Length  Proportional  to  Speed;  See  Maximum  Vector 
Displayed  in  I.ower  Right),  Terrain  elevations  are  highlighted  as  in  Figure  3 


Figure  8b.  Model  Windfield  for  Kxperiment  2  Showing  Wind  Velocity  Vectors  for 
Stable  Case  With  Vegetation  (Length  Proportional  to  Speed;  See  Maximum  Vector 
Displayed  in  Lower  Right).  Terrain  elevations  are  highlighted  as  in  Figure  3 
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Figure  8c.  Model  Windfield  for  Fxperiment  7  Showing  Wind  Velocity  Vectors  for 
Stable  Case  W  ith  Vegetation  (Length  Proportional  to  Speed;  See  Maximum  Vector 
Displayed  in  Lower  Right).  Terrain  elevations  are  highlighted  as  in  Figure  3 


3.3.2  UNSTABLE  CONDITIONS  WITH  VEGETATION 


The  experiments  simulating  unstable  conditions  were  experiments  3,  4,  5, 
and  8.  In  experiments  3  and  4,  the  input  wind  speed  of  2  ms  1  corresponds  closely 
to  the  climatological  value  (3  ms  *).  Use  of  the  same  wind  speeds  (2  ms  *  in  ex¬ 
periments  3  and  4)  enables  us  to  make  comparisons  of  these  experiments  with 
those  for  the  stable  cases. 

The  windfield  for  the  slightly  unstable  case  (experiment  3)  is  shown  in  Figure 
9a.  Notice  that,  for  this  case,  the  model  produces  upslope  flows  and  ridgetop  con¬ 
vergence  over  areas  of  high  terrain.  When  the  magnitude  of  the  buoyancy  is  in¬ 
creased  (experiment  4;  Figure  9b),  the  upslope  flows  are  enhanced,  and  some  of 
the  convergence  patterns  along  the  higher  terrain  are  sharper.  As  with  the  stable 
case,  the  model  is  capable  of  simulating  topographically  induced  flows  given  ini¬ 
tially  calm  winds.  This  is  illustrated  by  the  windfield  of  experiment  8  (Figure  9c). 
As  in  experiment  7,  the  strongest  slope  flows  nearly  coincide  with  the  steepest 
terrain  slopes. 

All  of  the  experiments  discussed  so  far  have  been  run  using  small  initial  wind- 
speeds.  With  such  light  winds,  it  would  be  expected  that  topographically  induced 
flows  would  be  the  dominant  features.  The  question  then  arises.  What  if  higher 
winds  prevailed  over  the  area?  One  would  expect  the  stronger  winds  to  dominate 
the  topographic  flows.  This  is  demonstrated  in  experiment  5  (Figure  9d),  which 
shows  a  fairly  uniform  windfield  only  slightly  influenced  by  topography.  An  un¬ 
stable,  high  wind  speed  case  can  exist  for  situations  where  vertical  mixing  of 
strong  momentum  from  aloft  begins  to  occur  after  breakdown  of  the  nocturnal  inver¬ 
sion.  It  has  been  shown  over  mountainous  regions  that,  once  this  strong  mixing 

has  begun,  topographically  induced  flows  begin  to  disappear  and  are  replaced  by 

30  3 1 

strong  uniform  wind  conditions  at  the  surface.  ’ 

3.3.3  NEUTRAL  CONDITIONS  WITH  VEGETATION 

Neutral  conditions  were  simulated  using  the  July  climatological  sounding  (Fig¬ 
ure  7c;  Table  5).  The  windfield  for  this  case  is  shown  in  Figure  10.  The  results 
are  as  expected:  fairly  uniform  winds  with  only  slight  influence  from  the  topo¬ 
graphy.  To  verify  that  the  uniform  windfield  was  caused  only  by  the  neutral  sur¬ 
face  layer  and  not  by  the  initial  wind  speed  of  10  ms"1,  an  additional  experiment 
was  run  with  neutral  conditions  and  an  initial  wind  speed  of  2  ms"1.  Results  of 
this  experiment  (not  shown)  revealed  that  the  uniform  windfield  was  indeed  caused 


30.  Banta,  R.M.  (1982)  An  Observational  and  Numerical  Study  of  Mountain 

Boundary-Layer  Flow.  Ph.  D.  thesis,  Colorado  State  U. ,  Ft.  Collins,  Colo. 

31.  Banta,  R.M.  (1984)  Daytime  boundary-layer  evolution  over  mountainous  ter¬ 

rain.  Part  I:  Observations  of  the  dry  circulations.  Mon.  Wea.  Rev. 
112:340-356. 
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Figure  9a.  Model  Windfield  for  F.xperiment  3  Showing  Wind  Velocity  Vectors  for 
I  nstable  Case  With  Vegetation  (l.ength  Proportional  to  Speed;  See  Maximum  Vector 
Displayed  in  Lower  Right).  Terrain  elevations  are  highlighted  as  in  Figure  3 


,  i 


28 


■ 

60- 70m  |  |  90 -100m 

70 -80m  El 

■  100 -110m 

FT  POLK 

TERRAIN  W/VEG  INCL 

■ 

80 -90m  I  |  >  110m 

20 

30 

40 

50 

i 


Figure  9h.  Model  Windfield  for  Fxperiment  4  Showing  Wind  Velocity  Vectors  for 
I  nstable  Case  With  Vegetation  (Length  Proportional  to  Speed;  See  Maximum  Vector 
Displayed  in  Lower  Right).  Terrain  elevations  are  highlighted  as  in  Figure  3 
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Figure  9c.  Model  Windfield  for  Experiment  8  Showing  Wind  Velocity  Vectors  for 
Unstable  Case  With  Vegetation  (Length  Proportional  to  Speed;  See  Maximum  Vector 
Displayed  in  Lower  Right).  Terrain  elevations  are  highlighted  as  in  Figure  3 


Figure  10.  Model  Windfield  for  Experiment  6  Showing  Wind  Velocity  Vector  for 
Neutral  Case  (Length  Proportional  to  Speed;  See  Maximum  Vector  Displayed  in 
Lower  Right).  Terrain  elevations  are  highlighted  as  in  Figure  3 


by  the  neutral  conditions.  Notice  from  Figure  10  that  speed-up  of  the  wind  occurs 

downwind  of  some  of  the  low  spots  in  the  terrain.  Speed-up  has  been  theoretically 
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related  to  terrain  geometry  (see  Orgill  for  a  short  review  of  the  work  done  on 
the  speed-up  factor),  so  it  is  not  surprising  that  the  model  is  capable  of  reproduc¬ 
ing  this  phenomenon,  if  somewhat  crudely. 

3.3.4  STABLE  CONDITIONS  WITHOUT  VEGETATION 

Experiments  simulating  the  stable  surface  layer  without  vegetation  and  with 
uniform  surface  roughness  of  14  cm  are  experiments  la,  2a,  and  7a.  Results  from 
experiment  la  are  shown  in  Figure  11a.  Downslope  flow  corresponding  to  cold  air 
drainage  appears  in  this  experiment.  Notice  the  uniformity  in  the  drainage  flow 
along  the  NW-SE  axis  of  the  Sixmile  Creek  valley;  the  drainage  flow  appears  to 
describe  a  "downslope  front"  coming  off  the  ridge.  Contrasting  this  flow  field 
with  the  one  in  experiment  1  (same  conditions,  except  with  vegetation  and  variable 
surface  roughness;  see  Figure  8a)  reveals  that  the  addition  of  vegetation  disrupts 
the  uniformity  of  the  downslope  front,  concentrating  the  drainage  instead  into  lo¬ 
calized  low  spots  corresponding  to  the  lowland  marshes.  These  results  have  sig¬ 
nificant  implications  in  the  area  of  toxic  chemical  diffusion.  The  presence  of  non- 
homogeneous  vegetation  can  alter  the  flow  field  and  subsequently  alter  the  chemical 
transport  and  diffusion,  perhaps  focusing  high  concentration  levels  into  these  lo¬ 
calized  lowland  areas  under  stable  conditions. 

When  the  buoyancy  is  increased  (experiment  2a),  similar  effects  take  place  as 
occurred  in  the  stable  cases  with  vegetation  (experiment  2).  The  windfield  for  ex¬ 
periment  2a  is  omitted  here  for  that  reason.  An  experiment  simulating  stable  con¬ 
ditions  with  initially  calm  winds  (experiment  7a)  was  conducted,  and  the  windfield 
is  shown  in  Figure  lib.  Notice  that  the  downslope  flow  appears  over  different 
parts  of  the  domain  than  it  does  in  the  same  case  with  vegetation  (experiment  7; 
Figure  8c).  This  is  a  further  illustration  of  the  effects  of  nonhomogeneous  vege¬ 
tation  on  drainage  flows,  again  with  important  implications  on  the  transport  and 
diffusion  problem. 

3.3.5  UNSTABLE  CONDITIONS  WITHOUT  VEGETATION 

Experiments  3a,  4a,  5a,  and  8a  were  simulations  using  unstable  conditions 
with  no  vegetation,  and  with  uniform  roughness  of  14cm.  The  windfield  for  exper¬ 
iment  3a  appears  in  Figure  12a.  Notice  the  upslope  flows  and  ridgetop  conver¬ 
gence  simulated  by  the  model  in  this  experiment,  consistent  with  the  results  of 
the  same  case  with  vegetation  (experiment  3;  Figure  9a).  However,  the  windfield 


32.  Orgill,  M.  M.  (1981)  Atmospheric  Studies  in  Complex  Terrain:  A  Planning 
Guide  for  Future  Studies.  PNI.-3656,  ASCOTr80/4,  Pacific  Northwest 
Laboratory,  Richland,  Wash. 
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Figure  11a.  Model  Windfield  for  Experiment  la  Showing  Wind  Velocity  Vectors 
for  Cases  With  No  Vegetation  and  Uniform  Surface  Roughness  (Length  Proportional 
to  Speed;  See  Maximum  Vector  Displayed  in  Lower  Right).  Terrain  elevations  are 
highlighted  as  in  Figure  2 


from  experiment  3a  contains  a  well  defined  southeasterly  component  along  the  Six- 
mile  Creek  valley  which  did  not  appear  in  the  experiment  with  vegetation  (compare 
Figures  12a  and  9a).  The  presence  of  nonhomogeneous  vegetation  in  experiment  3 
seems  to  disrupt  the  southeasterly  flow  along  the  valley,  with  the  lowland  marshes 
distorting  it.  The  southeasterly  flow  in  experiment  3a  appears  to  be  partially 
caused  by  the  channeling  effects  of  the  parallel  ridge  and  highlands  in  the  northern 
and  eastern  areas  of  the  domain.  We  do  not  believe  this  is  a  valley  wind  because 
a  corresponding  northwesterly  mountain  wind  does  not  appear  in  the  stable  case. 

When  the  magnitude  of  the  stability  is  increased  (experiment  4a),  the  south¬ 
easterly  flow  is  maintained  and  the  windfield  remains  essentially  the  same  as  in 
experiment  3a.  Therefore,  we  do  not  show  the  windfield  for  this  experiment. 

The  appearance  of  the  southeasterly  flow  along  the  valley  in  the  experiments  with¬ 
out  vegetation  would  have  important  effects  on  transport  and  diffusion.  The  re¬ 
sults  suggest  that,  in  the  cases  without  vegetation,  chemicals  could  be  transported 
further  up  the  valley  than  in  the  vegetation  cases,  where  diffusion  could  take  place 
in  the  localized  lowland  marshes.  When  initially  calm  winds  are  introduced  (ex¬ 
periment  8a;  Figure  12b),  the  localized  upslope  flows  appear  over  different  areas 
than  they  do  in  the  vegetation  case  (experiment  8;  Figure  9c).  When  an  initial 
wind  speed  of  10  ms'1  is  introduced  (experiment  5a;  Figure  12c),  the  windfield  is 
nearly  uniform,  more  uniform,  in  fact,  than  in  the  vegetation  case  (experiment  5; 
Figure  9d).  This  probably  occurs  because  the  terrain  is  more  uniform  without  the 
vegetation. 

3.3.6  NEUTRAL  CONDITIONS  WITHOUT  VEGETATION 

A  curious  result  comes  from  this  experiment  (experiment  6a;  Figure  13). 

The  windfield  for  this  case  is  totally  uniform,  remaining  unchanged  from  initial 
conditions.  The  effect  of  terrain  uniformity  on  the  flow  was  discussed  at  the  end 
of  Section  3.  3.5.  To  determine  if  the  uniform  windfield  was  the  result  more  of 
terrain  and  stability  conditions  than  of  initial  wind  speeds,  experiments  were  run 
with  neutral  conditions  and  an  initial  wind  of  2  ms  1  (no  vegetation),  and  other  ex¬ 
periments  were  run  with  no  terrain  and  uniform  roughness  for  differing  wind- 
speeds.  The  results  of  these  experiments  led  us  to  conclude  that  the  removal  of 
the  vegetation  and  imposition  of  neutral  conditions  was  responsible  for  the  uniform 
windfield  in  experiment  6a. 


i.  COINC.USIONS 

The  first  phase  of  the  ASL  windflow  model  evaluation  has  been  completed. 

We  have  shown  that  the  model,  despite  the  simple  physics  involved,  is  quite  ver- 
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Figure  12a.  Model  Windfield  for  Experiment  3a  Showing  Wind  Velocity  Vectors 
for  Cases  With  No  Vegetation  and  Uniform  Surface  Roughness  (Length  Proportional 
to  Speed:  See  Maximum  Vector  Displayed  in  Lower  Right).  Terrain  elevations  are 
highlighted  as  in  Figure  2 
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Figure  12b.  Model  Windfield  for  Experiment  8a  Showing  Wind  Velocity  Vectors  for 
Cases  With  No  Vegetation  and  Uniform  Surface  Roughness  (Length  Proportional  to 
Speed;  See  Maximum  Vector  Displayed  in  Lower  Right).  Terrain  elevations  are 
highlighted  as  in  Figure  2 
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Figure  12c.  Model  Windfield  for  Experiment  5a  Showing  Wind  Velocity  Vectors  for 
Cases  With  No  Vegetation  and  Uniform  Surface  Roughness  (Length  Proportional  to 
Speed;  See  Maximum  Vector  Displayed  in  Lower  Right).  Terrain  elevations  are 
highlighted  as  in  Figure  2 


satile  and  capable  of  simulating  different  types  of  flows  under  varying  topographic 
and  stability  conditions.  This  section  summarizes  the  modeling  results,  dis¬ 
cusses  some  of  the  limitations  of  the  windflow  model,  and  describes  the  future  re¬ 
search  with  the  model  that  will  be  done  at  AFGL. 

1. 1  Summary  of  Results 

The  ASL  windflow  model  can  simulate  stability-influenced,  topographically 
driven  flows,  such  as  nocturnal  drainage  and  upslope  flow.  The  model  can  simu¬ 
late  these  flows  under  light  wind  or  calm  conditions,  and  it  produces  fairly  uni- 
f  rm  windfields  under  high  wind  conditions.  Under  neutral  conditions,  the  model 
produces  essentially  uniform  windfields.  The  model  shows  a  tendency  towards 
intensifying  the  slope  winds  as  the  magnitude  of  the  buoyancy  is  increased,  while 
weakening  the  winds  elsewhere. 

Examination  of  the  total  kinetic  energy  (KE)  over  the  domain  for  different  sta¬ 
bility  conditions  revealed  several  insights  consistent  with  the  model  physics.  The 
KE  is  sensitive  to  the  magnitude  of  the  buoyancy,  increasing  faster  for  cases  with 
larger  buoyancy  magnitudes.  However,  this  result  only  appears  for  cases  with 
light  wind  conditions.  Under  high  wind  conditions,  the  KE  decreases  through  the 
model  integration.  Under  neutral  conditions,  the  model  KE  decreases  during  the 
model  integrations,  regardless  of  the  initial  wind  conditions.  These  results  are 
consistent  with  the  adjustment  taking  place  through  the  integral  of  Eq.  (1),  and  are 
thus  dependent  upon  the  buoyancy,  initial  wind  speed,  and  topography  in  the  model. 
The  strongest  winds  generally  appear  over  areas  with  the  greatest  terrain  gradi¬ 
ents. 

The  results  of  the  model  experiments  also  illustrate  the  importance  of  includ¬ 
ing  the  effects  of  vegetation  as  they  relate  to  the  modification  of  terrain  geometry. 
The  addition  of  vegetation  in  the  model  alters  the  topographically  induced  flows, 
creating  localized  convergence  zones  in  low-lying  marshes  under  drainage  (stable) 
conditions,  and  disrupting  the  southeasterly  channeled  flow  along  the  Sixmile  Creek 
valley  in  the  unstable  cases.  The  changes  in  the  flow  patterns  by  addition  of  vege¬ 
tation  are  consistent  with  the  model  physics,  since  the  vegetation  alters  the  ter¬ 
rain  slopes.  It  was  determined  that  variable  surface  roughness  is  a  second-order 
effect  in  the  model,  and  that  the  changes  in  the  flow  patterns  are  caused  primarily 
by  the  presence  of  vegetation  as  used  to  modify  the  terrain  geometry. 

1.2  Limitations  of  the  Research 

The  goal  of  this  work  is  to  produce  an  operational  prediction  system  for  toxic 
chemical  transport  and  diffusion.  The  numerical  model  described  in  this  report 
seems  well  suited  for  adaptation  to  an  operational  system,  but  it  has  some  impor- 


tant  limitations  which  should  be  pointed  out. 

The  simple  physics  of  the  model  as  described  by  Eq.  (1)  leaves  out  some  im¬ 
portant  terms,  such  as  surface-layer  momentum  fluxes.  Even  though  these  fluxes 
can  be  assumed  to  be  nearly  constant  in  the  surface  layer,  there  is  no  coupling  of 
the  surface  layer  to  the  mixed  layer  aloft.  This  is  an  important  weakness  in  the 
unstable  cases,  where  increasing  the  instability  should  lead  to  a  stronger  mixing 
and  more  uniform  wind  flow,  but  this  does  not  appear  in  the  model. 

The  vertical  dimension  is  only  considered  in  the  power  law  formulation  using 
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the  Businger  et  al  wind  profiles.  In  the  case  of  Fort  Polk,  it  can  be  assumed 
that  the  buoyancy  field  is  fairly  uniform,  but,  over  steeper  terrain  slopes,  effects 
of  differential  surface  heating  become  more  important.  In  the  latter  case,  a  ther¬ 
modynamic  equation  or  a  surface  energy  budget  might  have  to  be  added  to  consider 
the  effects  of  differential  surface  heating  and  vertical  mixing  of  heat,  momentum, 
and  moisture  on  the  surface  layer. 

Another  useful  addition  to  the  model  would  be  the  effects  of  the  vegetation  can¬ 
opy  itself  on  the  wind  profiles  and  stability  of  the  surface  layer.  Modifications  to 
the  ASL  model  to  include  vegetation  canopy  effects  are  presently  being  made  at 
ASL.  * 

Obvious  limitations  exist  on  the  availability  of  meteorological  observations  on 
a  scale  such  as  the  one  we  are  using  in  this  model.  This  problem  is  related  both 
to  the  availability  of  input  observations  and  to  the  verification  of  the  model  wind- 
fields.  The  model  experiments  conducted  in  this  evaluation  were  not  verified  by 

observations  over  the  Fort  Polk  area.  Some  data  were  collected  over  the  Fort 

19 

Polk  region  during  a  field  experiment  conducted  in  1980,  but  these  were  not  ade¬ 
quate  to  use  in  a  model  verification  study.  However,  the  results  presented  here 
are  reasonable  approximations  of  the  types  of  flows  one  would  expect  over  complex 
terrain  under  varying  stability  conditions. 

The  fact  that  the  terrain  only  varies  by  several  10s  of  m  over  Fort  Polk  does 

not  mean  that  topographic  flows  cannot  exist.  Observational  studies  by  Caughey  et 
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al  and  Lenschow  et  al,  and  modeling  studies  by  Brost  and  Wyngaard  show 

that  cold  air  drainage  and  slope  flo^s  can  occur  with  terrain  slopes  as  small  as 
33 

0.0014.  These  flows  were  also  shown  to  have  significant  effects  on  the  structure 
of  the  stable  boundary  layer.  The  fact  that  the  presence  of  tall  vegetation  can  alter 

*Cionco,  R.M.  Personal  communication. 

33.  Caughey,  S.J. ,  Wyngaard,  J.  C.,  and  Kaimal,  J.  C,  (1979)  Turbulence  in  the 

evolving  stable  boundary  layer,  J,  Atmos.  Sci.  36:1041-1052. 

34.  Lenschow,  D.  H.,  Stankov,  B.B.,  and  Mahrt,  L.  (1979)  The  rapid  morning 

boundary-layer  transition,  J.  Atmos.  Sci.  36:2108-2124. 

35.  Brost,  R.A.,  and  Wyngaard,  J.  C.  (1978)  A  model  study  of  the  stably  strati¬ 

fied  planetary  boundary  layer,  J.  Atmos.  Sci,  35:1427-1440. 
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the  terrain  slope  is  also  consistent  with  the  vegetation/no  vegetation  results  pre¬ 
sented  in  this  study.  The  model  results  concerning  the  effects  of  vegetation  on  the 
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flow  are  consistent  with  those  documented  by  Geiger  and  others.  Thus,  the  next 
step  in  the  model  evaluation  would  be  to  verify  the  windfields  against  actual  field 
observations. 

1.3  Future  Research  With  the  Model  at  AFGf. 

With  the  first  phase  of  the  model  evaluation  complete,  the  attention  at  A FG I . 
will  now  shift  to  verification  of  the  model  windfields  against  observations.  Future 
work  calls  for  running  the  ASL  model  using  simulated  and  actual  meteorological 
data  over  Vandenberg  AFB,  Calif.  Again,  detailed  terrain  information  (50-m  grid 
spacings  on  the  tebrain  data  tape)  is  available.  The  model  domain  will  be  chosen 
to  cover  an  area  the  size  of  the  southern  launch  complex,  where  meteorological 
tower  data  are  available  for  use  as  multiple  input  observations  and  for  verification. 
The  simulated  data  studies  will  use  the  same  format  as  the  Fort  Polk  study  (cli¬ 
matological  soundings)  in  order  to  study  the  sensitivity  of  the  model  to  stability 
and  the  steep  terrain  slopes  that  appear  over  the  Vandenberg  area.  The  tower 

data  experiments  will  likely  make  use  of  the  archived  data  base  available  for  this 
37 

area. 

The  proposed  work  intends  to  study  various  meteorological  situations  of  inter¬ 
est,  as  well  as  the  model's  sensitivity  to  the  number  and  spatial  distribution  of  in¬ 
put  observations  used.  The  results  of  the  proposed  work  will,  we  hope,  yield  re¬ 
sults  quantifying  the  real-time  input  data  requirements  for  the  model  and  reveal 
detailed  information  about  the  surface-layer  windflow  over  the  area  from  which  the 
Space  Transportation  System  (STS)  will  be  launched  in  the  future. 


36.  Geiger.  R.  (1965)  The  Climate  Near  the  Ground,  Harvard  University  Press. 

Cambridge,  Mass,  p.  356. 

37.  Schacher,  G.E.,  and  Stanton,  T.  P.  (1984)  Mean  Flow  and  Turbulence  in 

Complex  Terrain  NPS/ Vandenberg  Measurement  System.  NPS-61  -84-005, 


s-m# 


,v.  -- 


•  v 
•.*  < ■ 


0MNO 


'■/CO 

r 


References 


1.  Fleischer,  M.T.  (1980)  SPILLS  -  An  Evaporation/ Air  Dispersion  Model  for 

Chemical  Spills  on  Land,  Shell  Development  Company,  Westhollow 
Research  Center,  Houston,  Tex. 

2.  Kunkel,  B.A.  (1983)  Comparison  of  Evaporative  Source  Strength  Models  for 

Toxic  Chemical  Spills.  AFGL-TR -83 -0307,  ADA139431. 

3.  Kunkel,  B.A.  (1984)  An  Evaluation  of  the  Ocean  Breeze/Dry  Gulch  Dispersion 

Model,  AFGL-TR  -84-0313,  AD  A157165. 

4.  Ball,  J.A.,  and  Johnson,  S.A.  (1978)  Physically  Based  High  Resolution  Sur¬ 

face  Wind  and  Temperature  Analysis  for  EPAMS,  ASL-CR -78-0043-1.  ~~ 

(i.S.  Army  Atmospheric  Sciences  Laboratory,  White  Sands  Missile  Range, 
N.  Mex. 

5.  Holdeman,  J.T.  (1982)  Modeling  the  tropospheric  wind  field  over  complex 

terrain.  Symposium  on  the  Composition  of  the  Non  Urban  Troposphere, 
Williamsburg,  Va. 

6.  Dickerson,  M.H.  (1978)  MASCON  -  A  mass-consistent  atmospheric  flux 

model  for  regions  with  complex  terrain,  J,  Appl.  Meteorol.  17:241-253. 

7.  Dempsey,  D.  P. ,  and  Mass,  C.F.  (1983)  A  one-level  mesoscale  model  for 

complex  terrain,  Proc.  Sixth  Conf.  on  Numerical  Weather  Prediction,  Am. 
Meteorol.  Soc.  ,  Boston,  Mass.  343-347. 

8.  Haltiner,  G.J.  (1971)  Numerical  Weather  Prediction.  Wiley,  New  York. 

9.  Walmsley,  J.  L.  (1983)  Modeling  of  Boundary-Layer  Wind  Flow  Above  Com¬ 

plex  Terrain:  A  Review,  Report  ARQB-83  005-L,,  Atmospheric  Environ¬ 
ment  Service,  DownsvTew,  Ontario. 

10.  Lavoie,  R.L.  (1972)  A  mesoscale  model  of  lake-effect  storms,  J.  Atmos. 

Sci.  29:1025-1040. 

11.  Lavoie,  R.L.  (1974)  A  numerical  model  of  trade  wind  weather  on  Oahu,  Mon. 

Wea,  Rev,  102:630-637. 


12.  Fox,  D.G.,  Fosberg,  M.A.,  Marlatt,  W.  E. ,  and  Reeser,  W.  (1976)  Analysis 

of  mountain  air  quality,  Proc.  Third  Symposium  on  Atmospheric  Turbulence, 
Diffusion,  and  Air  Quality,  Am.  Meteorol.  Soc. ,  Boston,  Mass.  470-475. 

13.  Ohmstede,  W.  D.  (1981)  The  parameterization  of  battlefield  dispersion-new 

frontiers.  Workshop  on  the  Parameterization  of  Mixed -Layer  Diffusion, 

Las  Cruces,  N.  Mex.  279-287. 

14.  Ohmstede,  W.D.,  and  Stenmark,  E.B.  (1981)  A  model  for  characterizing 

transport  and  diffusion  of  air  pollution  in  the  battlefield  environment.  Work¬ 
shop  on  the  Parameterization  of  Mixed-Layer  Diffusion,  Las  Cruces, 

N.  Mex.  416-423. 

15.  Cionco.R.M.  (1981)  A  meteorological  approach  to  chemical  defense  over 

complex  terrain  with  vegetation.  Workshop  on  the  Parameterization  of 
Mixed-Layer  Diffusion,  Las  Cruces,  N.  Mex.  323-328. 

16.  Amlicke,  B.B. ,  and  Coleman,  I.  W.  (1984)  High  Resolution  Wind  (HRW)  Model. 

MRC/WDC-R-089,  Mission  Research  Corporation,  Alexandria,  Va. 

17.  I.anczos,  C.  (1970)  The  Variational  Principles  of  Mechanics  (4th  ed),  U.  of 

Toronto  Press,  Toronto,  pp.  106-110. 

18.  Dutton,  J.  A.  (1976)  The  Ceaseless  Wind,  McGraw-Hill,  New  York,  pp.  435- 

439. 

19.  Sadeh,  W.Z.,  Law,  F.W.,  Marlatt,  W.E.,  Fox,  D.G.,  Dietrich,  D.  L. ,  and 

Cionco,  R.M.  (1982)  A  Survey  of  Micrometeorological  Parameters  Within  a 
Forest  Canopy  at  Fort  Polk,  Louisiana,  ASL-CR-0100-1,  U.S.  Army  At¬ 
mospheric  Sciences  Laboratory,  W'hite  Sands  Missile  Range,  N.  Mex. 

20.  Bradshaw,  P.  (1969)  The  analogy  between  streamline  curvature  and  buoyancy 

in  turbulent  shear  flow,  J.  Fluid  Mech.  36,  Part  1:177-191. 

21.  Atkinson,  B.  W.  (1981)  Mesoscale  Atmospheric  Circulations,  Academic 

Press,  London,  pp.  13-20. 

22.  Panofsky,  H.A.,  and  Dutton,  J.  A.  (1984)  Atmospheric  Turbulence:  Models 

and  Methods  for  Engineering  Applications,  Wiley,  New  York,  pp.  124-125. 

23.  Monteith,  J.  L.  (1973)  Principles  of  Environmental  Physics,  Edward  Arnold 

Ltd. ,  London,  pp.  88-90.  — 

24.  Sellers,  W.D.  (1965)  Physical  Climatology,  U.  of  Chicago  Press,  Chicago 

pp.  150-151 

25.  Lettau,  H.  (1969)  Note  on  aerodynamic  roughness-parameter  estimation  on 

the  bases  of  roughness -element  description,  J,  Appl.  Meteorol.  JL828-832. 

26.  Tanner,  C'.B.,  and  Pelton,  W.  L.  (1960)  Potential  evapotranspiration  esti¬ 

mates  by  the  approximate  energy  balance  method  of  Penman,  J,  Geophys. 
Res.  65:3391-3413.  ~ 

27.  Businger,  J.A.,  Wyngaard,  J.  C.  ,  Izumi,  Y.  ,  and  Bradley,  K.  F.  (1971) 

Flux-profile  relationships  in  the  atmospheric  boundary  layer,  J.  Atmos. 

Sci.  28:181-18 9. 

28.  McNider,  R.T.,  and  Pielke,  R.A.  (1 984)  Numerical  simulation  of  slope  and 

mountain  flows,  J.  Climate  and  Appl.  Meteorol.  23:1441-1453. 

29.  McNider,  R.T.  (1982)  A  note  on  velocity  fluctuations  in  drainage  flows, 

J.  Atmos.  ScL  39: 1658 - 1 660. 

30.  Banta,  R.M.  (1982)  An  Observational  and  Numerical  Study  of  Mountain  Bound¬ 

ary-Layer  Flow,  Ph.  D.  thesis,  Colorado  State  I’.  ,  Ft.  CoTlins,  Colo. 

31.  Banta,  R.M.  (1984)  Daytime  boundary -layer  evolution  over  mountainous  ter¬ 

rain.  Part  I:  Observations  of  the  dry  circulations,  Mon.  Wea,  Rev. 
112:340-356. 


32.  Orgill,  M.  M.  (1981)  Atmospheric  Studies  in  Complex  Terrain:  A  Planning 

Guide  for  Future  Studies.  PNL-3656,  ASCOT/80/4,  Pacific  Northwest 
Laboratory,  Richland,  Wash. 

33.  Caughey,  S.J. ,  Wyngaard,  J.  C. ,  and  Kaimai,  J.C.  (1979)  Turbulence  in  the 

evolving  stable  boundary  layer,  J.  Atmos.  Sci.  36:1041-1052. 

34.  Lenschow,  D.H.,  Stankov,  B.  B.  ,  and  Mahrt,  L.  (1979)  The  rapid  morning 

boundary-layer  transition,  J.  Atmos.  Sci.  36:2108-2124. 

35.  Brost,  R.A.,  and  Wyngaard,  J.C.  (1978)  A  model  study  of  the  stably  strati¬ 

fied  planetary  boundary  layer,  J.  Atmos.  Sci.  35:1427-1440. 

36.  Geiger,  R.  (1965)  The  Climate  Near  the  Ground.  Harvard  University  Press, 

Cambridge,  Mass.  p.  356. 

37.  Schacher,  G.E.,  and  Stanton,  T.  P.  (1984)  Mean  Flow  and  Turbulence  in 

Complex  Terrain  NPS/ Vandenberg  Measurement  System,  NPS-61-84-005, 
Naval  Postgraduate  School,  Monterey,  Calif. 


